arXiv: 1509.02477vl [cond-mat.mtrl-sci] 8 Sep 2015 


Observation of the two way shape memory effect in an atomistic model of martensitic 

transformation 

E. A. Jagla 

Centro Atomico Bariloche and Instituto Balseiro (UNCu), 

Comision Nacional de Energia Atomica, (8400) Bariloche, Argentina 

We study a system of classical particles in two dimensions interacting through an isotropic pair 
potential that displays a martensitic phase transition between a triangular and a rhomboidal struc¬ 
ture upon the change of a single parameter. Previously it was shown that this potential is able 
to reproduce the shape memory effect and super-elasticity, among other well known features of the 
phenomenology of martensites. Here we extend those previous studies and describe the development 
of the more subtle two-way shape memory effect. We show that in a poly-crystalline sample, the 
effect can be associated to the existence of retained martensite within the austenite phase. We also 
study the case of a single crystal sample where the effect is associated to particular orientations of 
the dislocations, either induced by training or by an ad hoc construction of a starting sample. 

PACS numbers: 


I. INTRODUCTION 


Martensitic transformations are temperature driven 
non-diffusive phase transitions in which the atoms of a 
solid perform small individual atomic displacements be¬ 
tween a high temperature high symmetry austenite (A) 
phase, and a low temperature lower symmetry marten¬ 
site (M) phase [TT]. When the symmetry of the marten¬ 
site phase is a sub-group of that of the austenite phase, 
the system displays the remarkable shape memory effect 
(SME) [TJ [5] : The sample in the martensite phase can be 
mechanically deformed a big amount, yet, when the tem¬ 
perature is raised and the austenite becomes stable, the 
sample returns to its original shape. In simple terms the 
origin of the shape memory effect is associated to the fact 
that in the martensite phase, and even if the sample is 
mechanically deformed, each atom preserves the atomic 
neighborhood it had in the austenite, and when temper¬ 
ature is raised and austenite becomes stable again, each 
single atom returns to its original position. 

It is quite remarkable that in many cases, after a sam¬ 
ple was submitted to one (or more) of these cycles, it is 
able to “remember” the kind of deformation it was sub¬ 
mitted to, and upon a new cooling step, it deforms spon¬ 
taneously to the remembered shape, without the applica¬ 
tion of an external deformation. The shape change under 
this cooling-heating cycling is termed the two way (or all¬ 
round) shape memory effect (TWSME). The origin of the 
TWSME is more subtle than the simpler SME. It does 
not originate in general symmetry considerations, but 
in the remaining, after the initial cooling-deformation- 
heating training cycles, of atomic arrangements within 
the austenite that encode the form in which the sam¬ 
ple was stressed, and upon a new cooling, they favor the 
nucleation of martensite in preferential orientations. It 
has been observed that the TWSME can be triggered by 
small remaining pieces of the martensite phase within 
the austenite (this is referred to as retained marten¬ 
site), collections of dislocations with particular orienta¬ 


tions, or regions of the austenite with anisotropic rem¬ 
nant strains. [3ii[ia[iaii3iii] The generic scheme giv¬ 
ing rise to the TWSME is qualitatively depicted in Fig. 
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FIG. 1: Sketch of the shape changes in a sample display¬ 
ing a martensitic transformation. Martensite is schematically 
indicated by the dashed regions, which can be vertical or 
horizontal, indicating a martensitic variant elongated in the 
corresponding direction, (a) Starting anstenite sample, (b) 
Isotropic martensite after cooling, (c) Martensite deforms 
due to interconversion of variants, after an applied stretching 
along X. (d) Sample returns to the austenite phase after heat¬ 
ing, recovering totally (dl) or partially (d2) its original form 
(the amount or remaining deformation in (d2) may be hardly 
visible experimentally). This is the SME. In the last case, 
the retained martensite in (d2) can trigger the spontaneous 
shape change of the sample after cooling (e2), giving rise to 
the TWSME. 


Although well understood in qualitative terms, many 
details of the TWSME have not been totally elucidated. 
Particularly, the thermodynamical, or dynamical reasons 
that favor the formation of retained martensite or other 
kind of elastic distorsions that trigger the effect are only 
poorly understood. Also, other interesting questions 
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can be stated, as for instance to what extent a start¬ 
ing austenite sample can be engineered, in such a way 
to obtain a predefined deformation effect upon cooling, 
without the necessity of training cycles. 

In this paper we contribute to answer these questions. 
In section II we present an atomistic two-dimensional 
model that displays a martensitic phase transition be¬ 
tween a triangular and a rhomboidal phase. In section 
III we construct a polycrystalline sample with periodic 
boundary conditions, and apply a cooling-deformation- 
heating training cycle. After that, we investigate if the 
TWSME appears. We observe the appearance of the 
TWSME for the case in which the deformation during 
the training cycle is large enough, which is correlated 
with only a partial recovery of shape after heating in the 
training cycle, and with the existence of preferential ori¬ 
entations of retained martensite in the austenite phase. 
In section IV we show the same effect for the more re¬ 
alistic situation of a long bar upon bending. In Section 
V we study the TWSME in a single crystal sample with 
dislocations, showing how preferential orientations of the 
dislocations can justify the appearance of the TWSME. 
We also address the potentially promising possibility of 
designing a sample with an ad hoc distribution of defects 
that displays the TWSME. We show that a particular 
spatial distribution of dislocations can be proposed for 
which a robust and persistent TWSME is obtained. Fi¬ 
nally, in Section VI we summarize and conclude. 


II. MODEL 

Martensitic transformations typically occur due to 
an entropic stabilization mechanism of the austenite 
phasejS]. The bare interaction potential between the par¬ 
ticles has the deepest minimum at the structure of the 
martensite phase, whereas the austenite phase is at most 
a secondary metastable minimum. Yet, the basin of at¬ 
traction is wider for the austenite than for the martensite, 
so that thermal fluctuations generate a decrease of the 
free energy of the austenite with respect to the marten¬ 
site as temperature is raised. At a certain temperature 
the two structures can exchange stability, and the tran¬ 
sition occurs. 

A full numerical simulation of a martensitic transfor¬ 
mation requires the simulation of the dynamics during a 
long time, to observe the transition between the struc¬ 
tures taking place. In addition, the elementary time step 
of the simulation must be small enough in order to re¬ 
produce the thermal movement of the particles. On the 
whole, this fully first principle simulation requires a large 
temporal span that makes the simulation very time con¬ 
suming. Simulations along these lines have in fact being 
done[SH5], however the subtle effects we want to detect 
are hard to observe with this technique. 

In an alternative approach, we consider an effective 
interaction potential that already includes the effects of 
entropy. In this sense, this atomic potential must be ac¬ 


tually considered to be a free energy functional in the 
spirit of Ginzburg-Landau free energies used in general 
descriptions of phase transitions. We will assume that the 
effect of temperature is to change the form of this effec¬ 
tive free energy potential. The stability range between 
martensite and austenite structures is directly encoded 
in the form of the potential, and then the simulation will 
display the transition without the need to truly simulate 
the atomic thermal vibrations. 

We numerically study the system by solving for the 
time dependence of the particle coordinates according to 
a standard Verlet scheme in the presence of a frictional 
term proportional to velocity. This term takes out of the 
system any kinetic energy generated during the transfor¬ 
mation, and effectively keeps the system at a local energy 
minimum. 

The potential we use and its basic features were pre¬ 
sented in [9]. For completeness we give the form of the 
potential here. The interparticle potential V{r) is com¬ 
posed of several parts: 

V{r) = Vo + Vi + V 2 + Vs (1) 

where 

Vo = Ao[l/ri2-2/r6 + l] i/ r < 1 (2) 

Vi = [(r-l)2(r+l-2c)V(c-l)Vl (3) 

if r < c 

V 2 = -A 2 [(r - d 2 - S 2 )^(r - d 2 + S 2 )^]/s 2 (4) 

if d 2 - S 2 < r < d 2 + S 2 

^3 = ^3[(^ ~ ~ S3)^(r — ds -I- Ssj^J/sg (5) 

if ds - S 3 < r < ds + S 3 

and 0 otherwise in all cases. Vq is the repulsive part of a 
LJ potential and its weight in the total potential is mea¬ 
sured by the parameter Aq. The fourth order term pro¬ 
portional to Vi contributes with an attractive well to the 
total potential. The last two terms are fine tuning terms 
that provide a small minimum of amplitude A 2 centered 
at d 2 , and a small maximum of amplitude A 3 centered 
at ds- They were adjusted to penalize appropriately the 
triangular lattice, and/or favoring the martensitically re¬ 
lated structure. The potential is fully determined by 
the set of parameters P = {Aq, A2, A3, c, (i 2 ) S 2 ) dg, S3}. 
To study the triangular-rhomboidal transition (T-R) we 
use P = {Ao, 0.003,0.01,1.722,0.98,0.04,1.74,0.2} with 
variable Aq. Aq is the parameter that models the effect of 
temperature. There is a critical value of Aq ~ 0.067 such 
that for Aq < Aq (Aq > Ag) the martensite (austenite) 
phase is globally more stable. 

The lattice parameter of the austenite phase is approx¬ 
imately 0.995 at the transition point, increasing slightly 
with the value of Aq. The austenite phase is character¬ 
ized by two different interatomic distances that are ap¬ 
proximately 0.93 and 1.04 at the transition. In order to 
visualize parts of the sample in the austenite or marten¬ 
site phase, and since a direct examination of every atomic 
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position is impractical, we make the following analysis. 
Given a configuration of the particles, neighbor particles 
located at a distance compatible with the largest inter¬ 
atomic distance in the martensite phase are identified 
(for practical reasons we use a distance window between 
1.015 and 1.2), and those links are plotted as segments. 
This highlights regions of the sample in the martensitic 
phase. In addition, since our attention is focused on the 
shape of the sample, the following color code is imple¬ 
mented. When the angle between the long segment of 
the martensite and the x axis is less than 45 degrees, the 
segmant will be plotted in red, whereas if this angle is 
larger than 45 degress, it will be plotted as green. This 
color coding allows to identify at a glance if the sample 
(or a part of it) is preferentially stretched in the x or y 
direction. 


III. SHAPE MEMORY AND TWO-WAY SHAPE 
MEMORY IN A PERIODIC, POLYCRYSTALLINE 
SAMPLE 

The “shape” of a sample is a characteristic that is 
strongly connected to the existence of free surfaces. Yet 
as a hrst example, and since we are looking for tiny effects 
that must show up in the simulations, a different set up 
will be more convenient. We will study in this section a 
system with periodic boundary conditions within a rect¬ 
angular box of size lx, ly In this set up, the shape of the 
sample is characterized by the ratio Ix/ly The values of 
lx of ly are allowed to adjust dynamically during the sim¬ 
ulation according to the following mechanism. The inter¬ 
nal compressions fx and fy in the system in both direc¬ 
tion X and y are calculated along the simulation and are 
driven towards externally imposed values/“* and /“*, 
by the adjustment of lx of ly, according to the equations 


( 6 ) 

If fx^y are set to zero, the equilibrium values of lx,y de¬ 
termine the shape of a “free” sample. Instead, non-zero 
values of model the application of external forces 
along the x or y directions. 

The starting sample consists of a poly-crystal con¬ 
structed in the following way. A number of grain cen¬ 
ters and orientations are randomly chosen, and they are 
used to generate a given grain orientation around each 
of the centers. Each grain extends in the sample accord¬ 
ing to a Voronoi tessellation criterion, namely each sector 
of the sample is dominated by the nearest grain center. 
As a hrst step, this initial conhguration is relaxed with 
the true interatomic potential and Aq = 0.085, namely 
well within the stability range of the austenite phase. 
The values of lx,y are also allowed to relax during the 
process, under = 0. The conhguration obtained 
is shown in Fig. I^a). We observe that although the 
grains are in the austenite phase in bulk, in the grain 


boundaries there are atoms that are at the appropriate 
relative distance so as to trigger the martensite transfor¬ 
mation. As Aq is reduced and the martensite becomes 
progressively more stable, martensite crystals grow from 
the grain boundaries and invade the bulk of the grains, 
until all the sample is transformed [hg. [^b)]. In this 
transformation the form of the sample does not change, 
beyond some huctuations associated to hnite size effects. 



(a) (b) 


FIG. 2: (a) The starting polycrystalline sample after initial 
relaxation in the austenite phase. Plotting details and the 
meaning of the color code are provided in the text. Note 
that although the bulk of the grains are in the austenite 
phase, atoms with interatomic distances compatible with the 
martensite are widely present at the grain boundaries, (b) 
The sample after transformation to the martensite phase. 
There is no appreciable shape change during this transfor¬ 
mation. 

At this stage we apply an external deformation through 
the setting of non-zero values of /“* = —until we 
get a prescribed change in the shape of the sample, char¬ 
acterized by the ratio Ix/ly Fig. [^a) depicts the struc¬ 
ture of the sample after a 9% change in Ix/ly, and ulterior 
relieve of the external force and relaxation. Interconver¬ 
sion of variants caused by the external deformation (fa¬ 
voring x-elongated red ones) is apparent, and the sample 
retains a deformation of about 5%. From this configu¬ 
ration, the value of Aq is taken back to 0.085, and the 
final conhguration is shown in Fig. db). Note that al¬ 
though the sample mostly returns to its initial conhgu¬ 
ration, displaying SME, there are remains of martensitic 
variants, mostly red ones, that correspond to pieces of 
martensite compatible with the previously applied de¬ 
formation. This produces that the shape of the sample 
in Fig. [^b) is not exactly equal to that in Fig. Hb), 
but some elongation along x remains. Fig. [^c) and (d) 
show the equivalent situation that occurs when the ex¬ 
ternal stretching is performed along the y direction. In 
this case the remaining martensite in Fig. j^d) is mostly 
of the green type, and the sample retains an elongation 
along the y direction. 

The detailed evolution of the shape of the sample, as 
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FIG. 3: (a) The sample in Fig. ib) after external stretching 
in the x direction and further relaxation once the stretching is 
withdrawn, (b) State after re-transformation to the austen¬ 
ite. (c) and (d) The same for an initial compression along x 
direction. The black square is an outline of the sample in Fig. 
I^b), for comparison. Note that the remaining martensite in 
(b) and (d) corresponds to the deformation produced by the 
applied stretching. 


turn, this imbalance is responsible for the preferential 
growth of conveniently oriented variants when the sam¬ 
ple transforms to the martensite phase. 

It is worth to be mentioned that the TWSME we ob¬ 
serve is robust with respect to cycling between marten¬ 
site and austenite. Typically the effect is maximum in 
the first cycle, it diminishes about 20 % in the second cy¬ 
cle, and then it conserves this value for at least 20 cycles, 
that was the longest run we performed. 

In experimental realizations, it has been observed that 
the TWSME can be almost perfect: the sample may re¬ 
cover the deformed shape, even if the remnant deforma¬ 
tion is so small that it is undetectable. On the con¬ 
trary, the plot in Fig. |^c) shows that the amount of 
the TWSME in our case is tiny: it approximately corre¬ 
sponds to a spontaneous deformation upon transition to 
the martensite that is about one half of the remnant de¬ 
formation in the austenite after training. We think that 
this is related to the grain size of our sample. In fact, in 
additional simulations using smaller grains, the TWSME 
was almost unobservable. Instead, in a more demanding 
simulation using grains of linear size about three times 
those in Fig. the effect was much stronger, as the two 
crosses in Fig. |^c) indicate. The dotted line indicates 
in this case that the shape recovery is about five times 
larger than the remnant deformations. We think that 
this trend makes plausible that the TWSME we are ob¬ 
serving is compatible with that observed experimentally, 
where the grain sizes are much larger than those in our 
samples. 


measured by Ix/ly is presented in Fig. |^a): The shape 
does not change in the first A-^M transformation. The 
sample elongates upon the application of the external 
stress (S). When the stress is relieved (i?), the deforma¬ 
tion is reduced, but the sample remains deformed. When 
the sample returns to the austenite phase {A 2 ) it mostly 
recovers its original shape. This is the standard SME. 
From the configurations A 2 , we perform new austenite- 
martensite cycling, without applying any external stress. 
The evolution of the shape during this process is seen in 
Fig. {A 2 —>■ M 2 —t A 3 ...). We see that a systematic 
TWSME is observed, that follows the training imposed 
by the original external deformation. To better quantify 
this effect, we first plot in Fig. ib) the amount of rem¬ 
nant deformation after the first re-transformation to the 
austenite (A 2 ) as a function of the deformation caused 
by the external loading (i?). We see that this depen¬ 
dence is strongly non-linear. It is vanishingly small if 
the amount of applied deformation is low, but it rapidly 
increases at larger applied deformation. In (c) we see 
that the intensity of the TWSME is proportional to the 
remaining deformation after the training step indicating 
that preferential orientation of the remaining martensite 
within the austenite is the main responsible of TWSME. 
In fact, the remaining deformation in the austenite can 
be argued to be proportional to the imbalance between 
the different orientations of the retained martensite. In 


IV. SME AND TWSME IN A FREE BAR 

The results in the previous section show clearly that 
our model contains all the necessary ingredients to ex¬ 
plain the origin of the TWSME. It would be desirable 
however, to show the effect appearing in an experimen¬ 
tally achievable configuration, particularly in a sample 
with free surfaces. 

In order to do this, we consider the case of the bending 
of a bar shaped sample, with free surfaces. Bending is 
produced in the simulations by the application of external 
forces along the vertical (y) direction onto every particle, 
with a magnitude Fy that depends on the x coordinate of 
the particle as Fy = Fo((2a:/L)^ —1/3) (the bar is located 
between x = — L/2 and x = L/2), and Fq is chosen to 
achieve the desired degree of deformation. 

The internal structure of the sample is again poly¬ 
crystalline, but in this case the spatial distribution of the 
grains was chosen by hand. The reason for this choice 
is that randomly placed grains produce too strong vari¬ 
ations from sample to sample in our small systems. The 
outline of the grains, and the actual initial configuration 
after relaxation are shown in Fig. The system is sub¬ 
mitted to a sequence of transformations as displayed in 
Fig. | 6 l First: reduce Aq down to 0.057 to reach the 
martensitic state (a). In this process the shape of the 















5 



FIG. 4: (a) Evolution of the shape (Ix/ly) of the sample 
according to the following protocol. A: starting austenite 
(Fig. g a)). M: transformation to the martensite ([^b)). S: 
stretched sample upon application of different external loads. 
R-. relaxed sample after external force is set to zero (exam¬ 
ples in|^a),(c)). A 2 ’. transformation to the austenite, and 
observation of the SME. M 2 , S 2 and the following: shape 
change under cyclic austenite-martensite transformation, in 
the absence of external loading. This is the TWSME. (b) 
The remnant deformation as a function of the maximum de¬ 
formation. The form of this curve is clearly non-linear, (c) 
The amount of shape recovery in the TWSME as a function 
of the remnant deformation in the austenite phase. This re¬ 
lation is seen to be almost linear in the range of deformations 
analyzed, with a proportionality factor of about 0.5. The two 
crosses and the dotted line show a few results in a sample 
with larger grain sizes, where this factor is much bigger. 


sample does not change appreciably. Second: application 
of an external bending stress followed by the removal of 
the stress. The final configuration obtained is shown in 
(b), where it is seen that an important bending deforma¬ 
tion remains in the sample. Third: transform back to the 
austenite, by increasing Aq to 0.085 (c), in this stage the 
SME shows up, and the bar tends to recover its original 
straight shape. Note however that the recovery is not 
complete: a certain amount of deformation remains that 
is originated in pieces of a;-elongated (green) martensite 
in the lower part of the sample, and j/-elongated (red) 
martensite in the upper part. This imbalance of marten¬ 
sitic variants are the germ of the spontaneous bending of 
the sample when we transform to the martensite phase 
(Ao = 0.057) in the absence of external stresses. The 
final configuration is shown in (d). 



FIG. 5: The bar shaped sample. The black lines outline the 
grains in the starting sample. The actual configuration of 
the sample after initial relaxation in the austenite phase is 
shown (the crystalline orientation of each grain was chosen 
randomly). 



FIG. 6: Evolution of the bar along the treatment, (a) Free 
sample in the martensite state, (b) Bending caused by an 
external loading, after the loading has been relieved, (c) 
Transformation to austenite. The sample straights, show¬ 
ing SME, but an appreciable deformation remains, (d) Re¬ 
transformation to the martensite phase. The sample bends 
(a small amount) spontaneously. This is the TWSME. 


To quantify the TWSME, which is tiny for our small 
samples, we calculated the average vertical coordinate of 
the particles in slices of the bar at different horizontal 
coordinates, and these values are plotted in Fig. for 
the configurations in Fig. We see that in addition to 
the clear SME between (b) and (c), a small but clearly 
noticeable TWSME between (c) and (d) exists, in which 
the sample bends spontaneously without the application 
of any stress. This effect is persistent upon successive 
cycling between austenite and martensite, and it is of 
the same intensity as the one described in the previous 
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FIG. 7: Profiles of the samples in the previous figure, to quan¬ 
tify the degree of bending in each stage. 


section for a sample with periodic boundary conditions. 
In fact, quantifying the effect by the curvature of the 
bar, we see that the shape recovery when the sample 
transform to the martensite (the “TWSME” arrow in 
Fig.[7l) is about 0.4-0.5 of the remnant deformation (the 
(c) profile in the same figure), very much as in Fig. |^c). 


V. SME AND TWSME IN SINGLE CRYSTALS 
WITH DISLOCATIONS 

In the previous Sections, the origin of the TWSME was 
related to the remaining of martensite variants within 
the austenite structure. For the systems analyzed, these 
remains were favored by the poly-crystal nature of the 
sample, that allows for low energy atomic rearrangements 
at the grain boundaries. 

Experimentally, TWSME is known to occur also in the 
case of single crystal samples. In fact, although a per¬ 
fect austenite single crystal is not able to encode within 
its structure a memory of a previous deformation in the 
martensite phase, the unavoidable existence of defects 
-particularly dislocations- in the crystal structure pro¬ 
vides a mechanism to generate the TWSME. Dislocations 
are defects that produce anisotropic stresses in the crystal 
structure, which can favor the growth of martensite vari¬ 
ants of particular orientations upon cooling [iniiiiiis]. 

Experimentally, it is known that the TWSME can be 
induced in single crystals by training, this means that 
dislocations acquire a convenient distribution to trigger 
the shape transformation in some prescribed way, after 
the removal of the external stress. It is also possible to 
consider the case of a sample that is engineered with a 
collection of dislocations distributed in some appropriate 
manner, such that the sample may have a tendency to 
spontaneously elongate in one direction and contract in 
the perpendicular one upon the martensitic transforma¬ 
tion. Although at present to produce such an ad hoc 
distribution of dislocations is not feasible, this seems an 


interesting theoretical situation to consider. 



FIG. 8: Lattice structure around a single dislocation in a 
perfect triangular structure. The dislocation core is indicated 
by the ”T” symbol. Highlighted interatomic links indicate 
those that are in the range [1.015,1.2], well different than 
the equilibrium lattice parameter of the triangular structure 
(0.995). 

The elastic distortions that a single dislocation gener¬ 
ates on an isotropic elastic solid are well known. How¬ 
ever, it is not necessary to consider this deformation in 
full detail here. The following analysis is more appropri¬ 
ate to our purposes. We constructed by hand a single 
crystal structure with a dislocation in an otherwise per¬ 
fect single crystal sample. After preparation, the struc¬ 
ture was relaxed using the interatomic potential in the 
austenite phase (Aq = 0.085). The final structure is 
shown in Fig. As in previous figures, we also plot 
in that figure the interatomic links that are in the range 
appropriate to become the largest lattice parameter of 
the martensite structure, i.e., they will be the germs on 
which martensite crystals will grow. It is thus seen in 
the particular case of Fig. that this dislocation will 
induce a tendency to elongate the sample along the in¬ 
sertion line of the dislocation (y direction in this case) 
and a tendency to contract along the glide direction (x 
direction). 

In order to see to what extent dislocations are effec¬ 
tive to encode the TWSME, we first consider a simple 
case. To avoid the existence of free surfaces that can 
blur the effect of the dislocations, a single crystal with 
periodic boundary conditions is necessary. However, it 
is not possible to accommodate a single dislocation with 
these boundary conditions. Instead, we used a crystal 
with two dislocations as sketched in figure |^a). Note 
that both dislocations have the same orientation of the 
glide direction, thus contributing in the same way to a 
possible TWSME. The separation between the two glide 
lines was chosen by hand, however, note that the relative 
position of the two dislocation along this direction (the 
X separation of the dislocation cores) is adjusted by the 
system, as dislocations are very mobile along x. As the 
Aq parameter of the potential is reduced, the system re¬ 
mains within the austenite phase up to Aq ~ 0.060 (Fig. 
I^a)), well below the equilibrium transition value which is 
0.067. This indicates that there is an appreciable energy 
barrier to the nucleation of the martensite in the present 
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case. Reducing Aq further produces an abrupt transition 
to the martensite phase ([^b)-(e)), in which only two of 
the three possible martensitic variants are present. These 
variants are precisely those favored by the dislocations in 
the system, and generate a global contraction along x 
and expansion along y, as indicated in the last panel of 
Fig. [^by the outline of the initial and final shapes. The 
amount of the TWSME is in this case very close to the 
ideal value that can be expected in this case (about 7 % 
in It/I n) This is a clear evidence that the sample displays 
a strong TWSME. 






(e) 


(f) 


FIG. 9: Martensitic transition of a sample with two disloca¬ 
tions and periodic boundary conditions, (a) Stable configu¬ 
ration for To = 0.060. (b) to (e) spontaneous transformation 
for To = 0.059. (f) final configuration at To = 0.057. The red 
rectangle is the shape of the original sample in (a), and allows 
to see that there was an important shape change. Note that 
individual links between particles are not discernible in the 
main plots. The zoom in (f) shows a detail of the structure. 


The previous example shows -in a rather ideal 
situation- that the TWSME may be induced in our 
model by dislocations. Looking for the same effect in a 
more realistic set up, we tried to train a single crystal to 
display TWSME. Since some amount of disorder in the 
starting sample is necessary, we started with a sample 
with a small concentration of vacancies, and performed 
austenite-martensite cycles in the presence of an exter¬ 
nal stress. It was observed that the vacancies reacommo- 
dated in the form of dislocations, that were distributed in 
the sample in a more or less random fashion. However, we 


did not observe an over abundance of dislocations with 
the insertion line parallel to the stretching direction, that 
would favor the TWSME. In fact this was not observed 
beyond the limit of sample to sample variations. It is not 
clear to us why this attempt was unsuccessful. 


We thus decided to insert by hand an ad hoc distri¬ 
bution of dislocations that produce the TWSME in a 
sample with free boundaries. We concatenated many el¬ 
emental pieces similar to those in Eig. [^a), although the 
separation between dislocations was much shorter. The 
sample obtained, after a relaxation step in the austenite 
phase is shown in Fig. 10 a). A detail displaying the 
distribution of dislocations can be seen in Fig. lUta)- 
Note that all dislocations have the same glide directions, 
and all contribute to a shrinkage along x and expansion 
along y upon the martensitic transformation. The sam¬ 
ple after the transition to the martensitic phase is shown 
in Fig. 10’b), where an appreciable decrease of length 


is apparent. The sample returns to its initial length af¬ 
ter back-transforming to the austenite (Fig. [To|(c)). The 
shape change is persistent, as Fig. shows, although 
we also see in Fig. and that the sample de¬ 

grades after cycling. We suspect that the degradation 
effect would be lesser in larger samples, and also in a 
three dimensional set up where the movement of disloca¬ 
tions is more limited due to dislocation entanglement. 



FIG. 10: A free bar engineered with a collection of disloca¬ 
tions to induce the TWSME. (a) Structure after a relaxation 
with a parameter To corresponding to the austenite phase, 
(b) The sample after the first transition to the martensitic, 
and (c) back to the austenite, (d) The austenite sample after 
a few cycles, where some degradation is observed. Details of 
the configurations in the two outlined boxes are shown in the 
next figure. 


The present example is a theoretical verification that 
dislocations can be appropriately distributed in a single 
crystalline sample in order to induce the TWSME. 























FIG. 11: Detail of the configurations in the two square boxes 
highlighted in the previous figure. In (a), note the almost pe¬ 
riodic arrange of dislocations that were introduce by hand. In 
(b) we show a particular part of the sample that has suffered 
of an important degradation due to dislocation movement. 



FIG. 12: The shape of the free bar in Fig. as a func¬ 
tion of successive austenite-martensite transformations. The 
I parameter reported is proportional to the moment of inertia 
with respect to x, defined as 7 = {xi — Xi)^. 


VI. SUMMARY AND CONCLUSIONS 

In this paper we have used a tunable, classical, 
isotropic potential for identical particles in two dimen¬ 
sions to study the origin of the two way shape memory 
effect (TWSME) in martensitic transformations. Our 
model starts with an effective model at intermediate time 


scales, in which the actual atomic vibrations are not ex¬ 
plicitly resolved, but their effect is incorporated in an ef¬ 
fective change of the interatomic potential. In this way, 
subtle effects such as the TWSME can be studied with¬ 
out much numerical effort, something that would not be 
possible from a full ab initio simulation. 

We first studied the case of a polycrystalline sample, 
and made clear how the TWSME has its origin in the 
remaining of martensite pieces of the appropriate orien¬ 
tation inside the austenite, after a cooling-deformation- 
heating protocol. In addition, we were able to show that 
the same effect is observable in the more natural experi¬ 
mental configuration of a bar, under bending stresses. 

In a second part of the work we have seen how the 
presence of dislocations may induce the TWSME of sin¬ 
gle crystal samples, as dislocations produce a strain field 
on its neighborhood that favors the martensitic trans¬ 
formation into variants that produce a stretching of the 
sample in the insertion plane, and a contraction along 
the glide direction. We have seen that a properly chosen 
periodic arrange of dislocations is able to induce a bulk 
effect with an important shape change. 

These findings indicate, first of all, that the present 
model and simulation techniques are appropriate to dis¬ 
cuss subtle effects, as the TWSME, that are well un¬ 
derstood on a qualitative basis, but only rather super¬ 
ficially in their microscopical details. It also allows to 
study problems that are not immediately accessible to 
present experimental capabilities, as for instance to ad¬ 
dress the effect of a given dislocation distribution in the 
shape change of a given sample under martensitic trans¬ 
formation. All these result are promising, and give con¬ 
fidence that other even more demanding studies, as for 
instance the interplay of martensitic transformation and 
cracking and fracturing of a sample can also be studied 
with the kind of model studied here. 
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